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ABSTRACT 


Pulse responses have been used to calculate frequency response 
functions for several years. Questions have been raised whether use- 
ful pulse responses are obtainable from systems having backlash and 
saturation effects. It has been found in this investigation that by proper 
choice of pulse shape, duration and strength, it is possible to realize 


useful pulse responses despite system defects. 


First it is necessary that the pulse function have continuous deri- 
vatives to avoid infinite inputs. Secondly, the pulse strength must be 
adequate to proyide an output large compared to the backlash and other 
uncertainties. Finally the pulse should endure long enough to introduce 


a pulse of adequate strength without inducing saturation in the system. 


Pulse tests of two simple positional servo - mechanisms were 
performed. The AVSD system is characterized by angular velocity 
signal damping and the MAVSD system by modified angular velocity 


signal damping. Parameters pertinent to the investigation are as follows: 


(1) Backlash 0.25 to 1.0 milliradians 

(2) Gear Train 100 to 1 

(3) Motor FD 84-21-1 Diehl Manufacturing Company 
(a) maximum acceleration 9300 rad/sec 


(b) maximum velocity 415 rad/sec 
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(4) Pulse functions 

(a) Rectangular 

(b) Displaced cosine 
(5) Pulse strength 0.8 to 3.2 milliradian second 
(6) Pulse duration (PD)(NP)R from 0.1 to 1.0 


Rectangular pulse functions were found to be unsuitable because of the 
discontinuous nature of their derivatives. A displaced cosine pulse 
function of (PD)(NP)R = 1.0 anda pulse strength of three point two 
(3.2) milliradian seconds produced negligible response deviations for 


a backlash of one milliradian in both systems. 
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The object of this investigation is to dete rmine a rational basis 


for the choice of pulse functions to test systems with limited non-lin- 


earities. 
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CHAPTER 1 


INTRODUCTION 


Testing individual components and complete systems is a neces- 
sary phase of engineering. This practice not only enables verification 
of predicted performance, but also yields new information concerning 
properties of the system. Prediction of system performance relies on 
the existing knowledge of component and system properties. However, 
new inherent properties are frequently disclosed by the correlation of 
predicted and actual performance. Thus the complementary nature of 


these two aspects of testing is evident. 


Dynamic characteristics of systems and comporents are inves- 
tigated by analysis of their responses to specific input quantities. Of 
the many possible input quantities, two classes are most useful. These 
consist of transient inducing inputs and sinusoidal inputs. The relation- 
ship between the transient response to a pulse input and the steady- 
state sinusoidal response has been known for many years for systems 
associated with linear differential equations with constant coefficients. 
Practical methods for calculating the sinusoidal frequency response from 
the pulse function transient response were devised by Bromberg and 


others. 


Pulse function testing is particularly valuable when sinusoidal me- 
thods are impractical. For example, when the life of the system is short, 
as in guided missiles, sinusoidal responses cannot be measured di- 
rectly. Aircraft systems, being subject to changes in environmental 
conditions during test, do not submit readily to sinusoidal response 
measurement. However, the time required for pulse testing is much 
less than that required for sinusoidal techniques. Thus, pulse tests can 
be completed without materially affecting system life and before sig- 


nificant changes in environment occur. In other instances, damage to 
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equipment which can result from sinusoidal tests may be avoided with 
pulse function testing. Recommendations have been made in the interest 
of economy that pulse techniques be practiced in lieu of sinusoidal for 


even those occasions when sinuscidal tests are feasible. 


Objections have been raised to this recommendation. Statements 
have been made that physical systems having backlash in the gear 
train together with velocity and acceleration limi's of the motor wil. 
yield such poor pulse function responses that the method is useless. 
Simultaneously, it is asserted that this defect is not as important for 
sinusoidal tests. These objections were described by Mr. Rawiey 
McCoy of the Reeves Instrument Company, who suggested the subject 


of this thesis. 


From physical reasoning it would appear that pulse techniques 
should provide results at leas+ as good as those obtained from sinus- 
oidal tests in the presence of the non-linearities outlined above. These 
non-linearities exist in some measure in every physical System. 
Successful application of pulse techniques to physical systems are re- 
ported by Seamans, Blasingame and Clementson'*), Bollay!?), Gardner 
and Ross‘*) and others. When difficulties appear because of defects 
in the gear train or limitations of the motor, it seems reasonable to 
suppose that the difficulties may be overcome by the proper choice of 
the pulse function input. The purpose of this investigation is to furnish 


a rational basis for the selection of the pulse function. 


Physical systems are subject to interferences which prevent the 
controlled-member from performing inthe ideal manner. Backiash 
in the gear train is comparable to interferences in this respect. System 
design practices are to reduce the effect of interferences to a tolerable 
level described as the level of uncertainty. Similarly, backlash should 
be minimized, preferably to the level of uncertainty, before the system 
performs as desired. The strength of the pulse input must be great 
enough to make the uncertainty and backlash negligible quantities in the 
response. However, when the maximum value of the physical system 
response exceeds a critical value, saturation in the velocity or acceler- 


ation of the controlled-member occurs. The response magnitude is there- 
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fore limited if saturation effects are to be avoided. A range of response 
magnitudes is established by these considerations. between a lower limit 
derived from the backlash or uncertainty level and an upper limit imposed 
by saturation. The problem is to determine whether or not sufficient 
latitude exists in this range to provide necessary information for an 
analysis of the dynamic properties of a given s ystem. In this connection, 
the influence of pulse shape and pulse duration, as well as pulse strength, 
has been examined. A qualitative study of the distortion of the response 
resulting from velocity and acceleration saturation has been made to 


verify the existence of an upper limit on the response magnitude. 


A reasonable surmise is that the deviation from the ideal response 
arising from backlash will be no greater than that due to the general 
uncertainty level. The average deviation created by the backlash is ex- 
pected to be the smaller of the two. This is a consequence of the inter- 
mittant nature of backlash as opposed to the continuous nature of the level 
of uncertainty in a physical system. A consideration of this situation leads 
to the premise that the effect of backlash on the response will be indis- 


tinguishable from the uncertainty level. 


Pulse functions,in this investigation, have been applied to two 
selected systems to provide data for quantitative evaluation of the problems 
associated with pulse testing physical systems. These two systems 
have typical closed-chain structures. Selection of the systems was in- 
fluenced by the relative simplicity of their performance equations and 
the attendant numerical calculations involved. However, conclusions 


derived from these cases are believed to be quite general in character. 


The investigation reported in this thesis indicates that the res- 
ponses of physical systems derived from pulse tests will be as satisfac- 
tory as those obtained by sinusoidal methods. A definite basis for æ- 
lecting the proper pulse function is developed. A pulse shape such as a 
displaced cosine function is desirable. For this type of pulse function, 
it is found that pulse strength can be increased as the pulse duration 15 
increased. Thus it is possible to avoid saturation effects and yet res- 
trict the response deviation due to a limited amount of backlash. Under 
these conditions pulse testing of components and complete systems is 


certainly feasible. 
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CHAPTER 2 


DETERMINATION OF PERFORMANCE OF IDEAL LINEAR SYSTEMS 


An ideal linear system is one composed of a set of interacting 
variables so related that a variation of any one variable causes varia- 
tion of one or more of the other variables, and in which the principle 
of superposition applies. Linear systems may be ofthe open-chain, 


closed-chain, or multiple loop type. A simple closed-chain linear 


system is illustrated in Figure 2-1. 








Operating 
Component 











Feedback 
Component 





Fig. 2-1 Functional diagram of a closed-chain linear system 


The performance of an ideal linear system is defined by Draper, 


(5) 


tional relationships among the physical quantities present. These 


McKay, and Lees as the effects of the system in establishing func- 
physical quantities fall into two categories, namely system inputs and 
outputs. The inputs in general are independent of actions of the system, 
while the outputs are dependent both upon the inputs and the operation of 
the system. The performance of the system may then be described by 
the functional relationships between the output s and the inputs. The 
functional relationships are expressed mathematically for linear systenis 
in the form of kinematic and dynamic equations of motion of the system. 
These equations take the form of linear integro-differential equations 
with constanti coefficients. The coefficients are determined by the para- 


meters of the given system. 


i 


Draper et aN develop the concept of the performance function 
which may be derived from the integro-differential equation. The 
functional relationship between a given input and output may be repre- 


sented under zero initial conditions as shown below. 


q = [РЕЈ], , 9: 
out id [tin out] in 


In this expression qdo is the output quantity, q._ is the input quantity, 


БІЗ? is the ен performance function. The performance function 
embodies two components. The first, sensitivity, accounts for magnitude 
variations between inputs and outputs of the system. The second compcn- 
ent, called the response function, represents the time history of the 


magnitude variations between the inputs and the outputs. 


The form of the performance function of particular interesi is 
the steady-state sinusoidal relating function. This form describes the 
response of an ideal linear system to sinusoidal input variations. Sinus- 
oidal inputs may be represented by the rotating vector complex plane 
method or by a plot of the sine function with Wet as its argument and the 
amplitude of the variation as its coefficient. Both methods are illus- 
trated in Fig. 2-2. In this figure dinja s the amplitude, We is the an- 


gular forcing frequency, and t is time. 





Абија 





Фа Абпја 9 aud Qin* QCn)a sin “ft 


Fig. 2:2 Representation of sinusoidal input functions. 
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The output or response of the system will also be a sinusoidal function. 
The response will differ from the forcing function in amplitude and phase 
angle. Dynamic amplitude ratio (DAR) is the ratio of the response am- 
plitude to the input amplitude. Dynamic response angle (DRA) is the 
phase shift between the input forcing function (IFF) and the output 
response (OR). Dynamic amplitude ratio and dynamic response angle 
are both functions of the forcing frequency. Fig. 2-3 shows a typical 
steady-state sinusoidal response for an ideal linear system. The steady- 


state sinusoidal relating function may then be expressed as: 


аа 


[РЕ]; | БІТІМІН | 
id faini Tout] Id E Jout 


The steady-state sinusoidal response may be expressed as 


= “out 


jw; 
Ч(їп)а ы: 


dout > Im [id 


out 








RESPONSE 


Sout  З(ош)а sin] wet AOR [4 ) 401] 


Qin "Qm a Sin "£1 


(їп) а 


Fig. 2-3 Representation of the steady-state response of an ideal linear 


system to a sinusoidal input. 
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The variation of the steady-state relating function with the forcing 
frequency is the frequency response of the system. Nondimensional 
performance expressions are derived by dividing the dynamic ampli- 
tude ratio by an arbitrarily selected reference amplitude ratio. Per- 
formance function derivations and frequency responses for first and 
second order ideal linear systems are contained in Instrument Engineering”, 


A typical linear scale nondimensional frequency response is shown in 


Fig. 2-4. 


ЕС 


PHASE ANGLE Ó IN DEGREES 





NONDIMENSIONAL. AMPLITUDE RATIO 


FREQUENCY RATIO (FR)= 6 = при Пиве 


Figure 2-4 Linear scale sinusoidal response characteristics. 


Pulse functions may be used as forcing functions for ideal linear 
systems. Two types of pulse forcing functions are considered by the 


authors. 


2-1 Rectangular pulse function (RPF) 


Pulse functions of this type are constant for the duration of the 
pulse. Discontinuities exist at the beginning and at the end of the pulse. 
The pulse is described by three attributes. These are pulse height, 
(PH), pulse duration(PD, and pulse strength(PS).To mathematically 
define the rectangular pulse function, let u(t) represent the pulse as 


a function of time, t. 
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Е ЕСО 
u(t) - (PH); 0*t $ (PD) 
u(t) = 0 ; t» (PD) 


Fig. 2-5 shows a representation of a rectangular pulse function 


(PH) 


О PD ТІМ Е њи 


Fig. 2-5 Representation of a rectangular pulse function. 


Pulse strength (PS) is defined 


(PD) 
(PS) «f (PH)dt - (PH)(PD) 
4Q 


„ж 


2-2 Displaced cosine pulse function (DCPF) 


Pulse functions of this type vary in magnitude throughout the 
duration of the pulse. They attain a maximum value only once during 
the pulse duration. The pulse is described by three characteristics. 
These are pulse height (PH), pulse duration (PD), and pulse strength 
(PS). To mathematically define a displaced cosine pulse function, let 


u(t) represent the pulse as a function of time, t. 
ІШІ 9:9 (20 


u(t) cm. - cos mi 0 ту 
2 (PD) 


ut) 0; t > (PD) 
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Fig. 2-6 shows a representation of a displaced cosine pulse function. 


PH) 


Fig. 2-6 Representation of a displaced cosine pulse function. 


Pulse strength (PS) is defined: 


PD 
Es JA mu ы 3 at - (PHPD) 
о 2 (РГ) 2 


Instrument Engineering, ?) Bromberg Seamans, Blasingame, and 
Clementson Запа others develop numerical and graphical methods for 
obtaining the frequency response of an ideal linear system fro m data 
obtained through pulse testing. Essentially the methods consist of apply- 
ing a pulse disturbance to the system and recording the input pulse and 
output response as functions of time. Fig. 2-7 shows illustrative time 
responses of second order (0;0, 1, 2) linear systems for pulse function 
inputs. The required performance function form may be determined from 
the time histories recorded. The procedure is to obtain the Fourier 
transforms of the input and the response. These transforms may be 
shown to be continuous functions of frequency. The ratio of the Fourier 
transform of the response to that of the input is the desired relating 
function. The relating function thus determined expresses the complete 


frequency response of the system. 
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О PD + — 


(b) Response to displaced cosine pulse function 


Fig. 2-7 Linear second order system time responses to pulse function 
inputs 


From the foregoing it is evident that for an ideal linear system 
the same criteria may be determined by pulse disturbance methods 
as can be determined by sinusoidal methods. This is noteworthy since 
the time required to apply pulse test techniques is much less than that 
required using sinusoidal methods. Further it is conceivable that pulse 
techniques can be applied to systems when sinusoidal testing might cause 
excessive wear or damage to the system whose performance is to be 


determined. 
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CHAP T Hives 


DESCRIPTION OF THE PERFORMANCE OF SYSTEMS WITH LIMITED 
NON-LINEARITIES 


A non-linearity exists in a system when a dependent variable is 
determined by an independent variable in accordance with some non- 
ideal relationship. For example, the independent variable might vary 
sinusoidally and the dependent variable will then vary periodically but 
not sinusoidaily. This constitutes non-linear dependence. Fig. 3-1 


illustrates the non-linear relationship just described. 


делу > independent variable 
Giout) =dependent variable 





mae. 3-1 An example of a non-linear response. 


The non-linearities considered in this thesis are nominal amounts of 
backlash in a gear train as well as saturation characteristics ofa 
motor. As noted in Chapter 1, these will exist in some measure in 


almost every physical system. 


A system response forced by a pulse function input is termed a 
transient response. An example of such a response is shown in Fig, 2-7 
of Chapter 2. The transient nature is readily evident. Correspondence 
of the input and output quantities is attained only when both have returned 


to zero displacement. The response does not reach a state of continuous 
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correspondence with the input pulse function prior to this time. Further, 
it does not follow the pulse function with a constant error. The discussion 


which follows is concerned only with system transient characteristics. 


Fig. 3-2 illustrates schematically the relationship existing between 
the controlled-member and the torque summing member when backlash 
is present. The term controlled-member refers to the output member 
of a feedback control system. The torque summing member is considered 


to be the next to last gear in the train. 


Tustin ^) develops a technique that makes it possible to closely 
approximate the controlled-member motion when backlash exists in the 
gear train. It is assumed that some point between the torque summing 
member and the controlled-member moves in the manner predicted by 
ideal linear theory. This assumption is reasonable when the backlash 
constitutes a relatively small part of the total motion of the controlled- 
member. The location of the point is established by Tustin's method 
as the mean value of the motions of the controlled-member and the torque 
summing member, Motion of this point is referred to as mean motion. The 
actual motion of the controlled-member is calculated by combining the 
mean motion predicted by ideal linear theory with a deviation calculated 
from the backlash. The angular displacement of the point performing 


the mean motion is given by the expression: 


A Е Кет) А (ст) i ltsm) ^ (tsm) (3-1) 


(mean) 
(ста) ш l(tsm) 


where 


ltem) controlled-member inertia 


I - torque summing member inertia referred 
(tsm) 
to the control member. 
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Fig. 3-2. Schematic diagram of the relationship between controlled - 
member and torque summing member when backlash exists 
in the gear train. 
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The motions of the controlled-member and torque summing mem: 
ber are intermittant relative to the mean motion as the backlash is taken 
up in one direction or the other. When the backlash is small the motion 
of each member may be represented as a deviation from the mean motion. 
From the definition of the mean motion, it is noted that the deviation of 
the torque summing member is opposite to that of the controlled- member 
when the backlash is being taken up. The sum of these deviations equals 
the whole backlash. Each member may be considered to take up a fixed 
proportion of the total relative motion due to backlash. The portion of 


the backlash taken up by the controlled-member is: 


l(tsm) P (bl) 


= А 3-2 
ТН 2. (bl)(em) Ше) 
(tsm) (cm) 
That taken up by the torque summing member is: 
Қош) Мы) 2, 
(bl)(tsm) (3-3) 


ltem) Ы ст) 


It is now possible to describe the motion of the controlled member. 
Referring to Fig. 3-3a, the decaying transient oscillatory curve repre- 
sents the velocity of the mean point. The backlash is assumed to be 
taken up at the initial instant. As a result, the controlled-member vel- 
ocity duplicates that of the mean point until the first peak occurs. At 
this instant the mean velocity decreases. However, in the absence of 
friction, the controlled-member continues to move with a constant velocity 
because of the backlash. The mean velocity decreases until the effective 
backlash is taken up at which instant the gear teeth make contact. It 
is assumed that contact between the gear teeth is achieved without bounce 
effects. At contact the torque summing member slows down the controlled- 
member wnile it in turn speeds up. Fig. 3-3b shows the velocity devia- 


tions of the controlled- member relative to the velocity of the mean point. 
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Fig. 3-3a Velocity of the controlled-member relative to the velocity 
of the mean point for backiash in the gear irain. 
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Fig. 3-36. Deviation velocity of the controlled-member 
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Fig. 3-3c. Displacement of controlled-member relative to the 
displacement of the mean point, 


Figure 3-3. Motion of the controlled-member relative to the mean point 
for a decaying transient osciliarion,of the mean point as 
adapied from Tustin” 


oul 


The controlled-member displacement curve is obtained by integrating 
the velocity curve of the controlled- member. Fig. 3-3c represents 

the angular displacement curves corresponding to the velocities in Fig 3-3a. 
The integral of each of the shaded triangular areas of Fig. 3-3a is equal to 


the effective backlash experienced by the controlled- member. 


From the above discussion, it is apparent that the deviation of the 
controlled- member displacement never exceeds the effective backlash 
of the controlled-member. When the backlash itself is small and the 
peak displacement of the mean point is comparably large the net effect 
of backlash is insignificant. For this case, the area between the con- 
trolled-member and mean point angular displacement curves will be a 


very small percentage of the area under the mean point curve. 


Torque saturation characteristics of the motor may be considered to 
impose limits on the velocity and acceleration of the controlled-member. 
The undesirable effects of driving a system into saturation may be seen 
in Fig.3-4. This figure shows the ideal linear angular displacement 
response of a second order system to a rectangular pulse input. Res- 
ponse curves are also included for both velocity and acceleration 
limiting of the controlled-member. It is noted that velocity and accelera- 
tion limiting increase the peak amplitudes of the transient response. 
Further, a considerable time lag resulting from the limiting is evident. 
Fig. 3-5 illustrates the effect on the velocity of the controlled-member 
of velocity limiting and acceleration limiting. Similarly, Fig. 3-6 
shows the effect of acceleration limiting on the acceleration of the con- 
trolled-member. Certain discontinuities are pointed out in the curves 
of Figs. 3-5 and 3-6. These are attributed to the effect of the discon- 
tinuity in the rectangular pulse function which occurs when it returns to 
zero. These curves were obtained from the General Purpose Simulator 
of the Instrumentation Laboratory of M.I.T, 


Qualitative study of Fig. 3-4 indicates that distortion of the transient 
response resulting from either velocity or acceleration saturation is 


unacceptable. The frequency response calculated by Bromberg's (1) 
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method from such a transient response will also be distorted. Satis- 
factory correlation between the predicted and actual system response 
under these conditions would be difficult, if possible. It is concluded 
therefore that the pulse strength must be limited to avoid saturation 


effects for successful application of pulse test techniques. 


The character of the pulse shape, as well as the pulse strength, 
becomes important if the controlled- member velocity and acceleration 
are to be prevented from reaching saturation. In this connection, Deri- 
vation Summary 3-1 develops the initial accelerations required of repre- 
sentative second and third order systems for a step input. It is evident 
that a step input to the second order system will demand an acceleration 
of the controlled-member at the initial instant equal in magnitude to 
the square of the undamped natural frequency of the system multiplied 
by the magnitude of the step. For the third order system, the same input 
demands an acceleration approximately equal to one and three tenths 
times as great. The argument may be extended to show that a rectangular 
pulse will force such accelerations at both the beginning and end of the 
pulse. Consideration of these facts suggests that acceleration, rather 
than velocity limitations are critical for rectangular pulse functions. 


The problem clearly will not be alleviated by imreasing the pulse duration. 


A displaced cosine pulse function, in contrast to the rectangular, 
possesses no discontinuities. Therefore it is anticipated that the relative 
importance of the velocity and acceleration limits will depend both on pulse 


strength and pulse duration. 


The discussion presented here is concerned with a method for cal- 
culating the angular displacement response of a system with backlash in 
the gear train. It is suggested that the response deviations from ideal 
will not exceed the backlash provided velocity and acceleration saturation 
are not encountered. Acceleration saturation will be difficult to avoid for 
rectangular pulse inputs. For small amounts of backlash, it is presumed 
to be possible to obtain sufficiently large peak response amplitudes with 
a displaced cosine input while avoiding saturation. If this can be done, the 
non-linearity of the transient response of the physical system may be 
considered negligible. Close correlation of the predicted and actual system 


responses then becomes practicable. 
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Consider the following second order linear differential equation 





with constant coefficients. 


1 2 2 
І = ! 
(p^ 2$ w, p tw. оу ЛЛ (ан (1) 
мћеге 
'р = non-dimensional derivative operator 
А,. ү = independent variable 
(in) 
А tout) dependent variable 
Wa = non-dimensional undamped natural angular frequency 
УУ damping ratio 
Т  - dimensional similarity factor 
The Laplace transformation yields: 
o 
а WA 
А (оці) (А) А ап) 09 (2) 


(On? ar 2 w At Ww) 


гог А= + ] 218 
andAÀ- nondimensional Laplace transform variable 
X= real component of À 
27/5 = imaginary component of À 


D 





ee A, - unit step function 
(in) 
| 
Тћеп, А (51) CÀ) UU У (3) 
Substituting in (2) 9 
6 UN 
р 4 . 
ony ae 2 (4) 


t | 
МА + 2ў!м At Wa) 


Derivation Summary 3-1 Initial acceleration of second and third order 
systems for a step input. (Page 1 of 3) 
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у ом У ДЕ 
A = 
ПТ, A 
(оп) E А(А +2} А ) 
g^ у 2 ^ 


A" у) = 
(out) A) Tu 2$'w A+ м 2) 


For simplicity let Wa е] ОУ o = 1.0 


Substituting for 'W andg, and expanding in partial fractions 


AY un (A) EM. SE 
C ны  — A ee 
(out) 
2j i -37 (A-N) 2,1 -92 (2-6) 


Where Ая is the complex conjugate of А, . Taking the inverse transform, 


A" E 9+1 2 At А у= 1 42. Ж" 





(7) 


= е E 8 
(out) 7-9 П? (8) 
231 - $ Zin tS 
where 't = non-dimensional running variable. 
At 't = 0 
А" = 1.0 = 'w - 9 
(out) | n (9) 


(a) initial acceleration of a second order (0;0,1,2) system for a step input. 
The third order equation similar to (6) is as follows: 
+ 1 
AX uty) = E (10) 
(А +2Д+1)(1+нА) 


мћеп E = 077 1.0 


t1 and to are first order non-dimensional time constants. 





Derivation Summary 3-1 Initial acceleration of second and third 
order systems for a step input. (Page 2 of 3) 
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The inverse transform gives: 
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Тпеп, А" У їогЧ = 0 (12) 


= 2 
RC Aon 1-22 


(0) Initial acceleration of a third order (1,0;0,1,2,2) system fora 


— == 





step input. 





Derivation Summary 3-1 Initial acceleration of second and third 
order systems for a step input, (Page 3 of 3) 
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CHAPTER 4 
GENERAL DESCRIPTION OF TWO SELECTED PHYSICAL SYSTEMS 


In Chapter 3, a method is proposed for determining the perfor- 
mance Of systems with some non-linearities. In order to find out how 
practical the techniques are, two physical systems are selected for 
study. Each is a positional servomechanism with two feedback branches. 
The first is one associated with a second order (0;0,1,2) performance 
equation. The other is a system associated with a higher order perfor- 
mance equation, namely a third order (1,0;0,1,2,3) type. Both systems 
have as an input angular displacement - А (ап) The output of each sys~ 
tem is a controlled-member angular displacement "(em 


4-1 Ideal Linear Properties of the AVSD and MAVSD Systems 


Д— —— = 


The first system is an angular velocity signal damped positional 
servomechanism which will be referred to as the AVSD system. A 
functional diagram of the AVSD system is shown in Figure 4-1. Its 
ideal linear performance equation is a linear differential equation of 


the following derivative operator form: 


~ 
= 


cn 


2 + | E | 
\Щста)ев Р у (рэ)ГА; мр} P (em) j P(ps) A: M] M 


(С) А intfr 


(ст) ” 
(4-1) 
The derivation of the above equation and definitions of the coefficients 


may be found in Appendix A. 


The ASVD system is characterized by the means by which damping 
is achieved. Negative unity feedback of !he controlled- member angular 
displacement causes generation of a correction signal whenever the dis- 
placement is different from that of the input, The existence of signal 
level elements in the chain of Fig. 4-1 permits reduction of the correc- 


tion signal by a signal proportional to the torque summing member angular 
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апо) у | шә іше TED В (юр)ә [ор азы 56 (од) а E (ме) о Е | «Dv 


= —(PSAP)(Sd) WSINVHOAWOANZS TWNOILISOd G3dWVd лумена ALIDOTSA UvINONY— — — — — 1l 


velocity. The reduced correction signal produces an actuating torque 
in the torque generator which tends to drive the correction signal to 
zero. Under these conditions, the actuating torque is the diffevence 
between a torque proportional to the original correction signal and a 
torque proportional to the angular velocity of the torque summing 
member. Thus a damping-like term, introduced at the signal level, 
is included in the actuating torque. This technique exactly simulates 


viscous damping while avoiding high power dissipation. 


As a final comment on the AVSD system, it can be shown that an 
error will exist between the output and input shafts when the angular 


velocity of the controlled-member is constant. 


The second system studied is shown in the functional diagram of 
Fig. 4-2. The diagram differs from that of Fig. 4-1 by the addition 
of a velocity signal modifier in the angular velocity feedback path. 
Therefore this system will be referred to as the MAVSD system 


hereafter. 


The velocity signal modifier consists of a passive network differ- 
entiating circuit. Incorporating the velocity signal modifier gives the 
following derivative operator form of the system ideal linear performance 


equation. 


ж 


(СТ), sm P 


Lr (СТ) sm)? 


2 


Ж ЗИ? | 
(emea) * 89) А, КЕСЕДІ *(tb)(av)[A; M]P (епу 


M (4-2) 


З(рв)[А;м] ©) (сет) * Minter. 
This equation is equivalent to a third order linear differential equation 
with constant coefficients. The reader is referred to Appendix A for the 
derivation of the performance equation and the definition of the coeffi- 


cients. 


A frequency dependent damping characteristic results from the 
inclusion of the velocity signal modifier in the MAVSD system. The 
output of the angular velocity signal generator is supplied to the differ- 


entiating network of the velocity signal modifier. The resultant output 
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is then combined with the correction signal senerated from negative 


unity feedback of the output position. 


The frequency -dependence of the velocity signal modifier is appar- 
ent in the steady-state frequency response form of its performance func- 


tion. 


(CT) gm) wp 


EL (vsm) e;e 
TO M) (CT) as У 


With the characteristic time equal tothe undamped natural period, it 

is evident that the function approaches zero for small values of the 
forcing frequency. For values of the forcing frequency of the order of 
the undamped natural frequency, the function approaches unity. The 
undamped natural frequency referred to isthat of the second order system 


existing before the inclusion of the velocity signal modifier. 


The desirability of frequency-dependent damping is developed by 


(8) 


teristics should not be prominent for low frequencies, but should become 


Lees Briefly, it is shown in the reference cited that the damper charac- 
important for frequencies in the range of the undamped natural frequency. 
This requirement is met by the velocity signal modifier of the MAVSD 


V 
7 


system. 


= 


4-2 Physical Properties of the AVSD and MAVSD Systems 


So far no consideration has been given to the physical properties 
of these two systems. Generalities have been purposely used to ensure wide 
applicability of possible corclusions. It must be realized that equations 
4-1 and 4-2 are not associated solely with the systems shown in Figs 
4-1 and 4-2. A variety of components may be chosen. For example, 
a large range of sensitivities and gear ratios is available. Equations 
4-1 and 4-2 are still applicable to describe the performance of the systems. 


This demonstrates the power of this method of investigation. 
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To obtain quantitative results from which conclusions may be 
drawn, it is necessary to specify numerical values of the parameters 
involved. Appendix A lists the numerical coefficients of the equations 
of motion of the systems selected for study. These coefficients are 


derived in Appendix A from physical properties of the systems. 


The physical systems chosen use a motor to drive a controlled- 
member in response to an input command. To drive the controlled- 
member, a gear train is often needed. The presence of the gear train 
introduces backlash into the system. In addition interference torques 
or forces are present. The effects of these ph ysical factors must be 


considered. 


It is known that backlash in a system causes a steady state as well 
as a transient displacement error. It is apparent that in a well designed 
system this error must be held to a minimum. To describe this error it 
is necessary to establish a representative value of backlash to be found 
in physical systems. Such a value may be dete rmined by considering 
Fig. 3-2. It is seen that the amount of backlash is the gap between ad- 
jacent teeth when gears are meshed. Angular backlash may be expressed 
in terms of this gap divided by the radius of the gear. Current machining 
tolerance may be specified in ten-thousandths of an inch. Therefore 
the gap may be as large as four times the tolerance. Assuming a toler- 
ance of plus or minus one ten-thousandth of an inch and a gear radius 
of two inches, it is seen that a backlash of approximately one half a 
milliradian is representative. Therefore in this study the range of 


backlash is one-fourth to one milliradian. 


Motors have limitations on output torque, angular velocity, and 
angular acceleration. Each of these limitations must be considered. 
Limited torque specifies a lower limit on the gear ratio that may be 
used. If too low a gear ratio is specified, the mechanical advantage 
of the motor will be insufficient to move the heavier controlled-member. 
A value of gear ratio of one-hundred to one from the motor to the con- 
trolled- member is chosen as typical. The value of gear ratio selected 
limits the controlled -member velocity and acc eleration to one one- 


hundredth of the motor maximum velocity and acceleration if saturation 
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is to be avoided. Chapter 3 explains why saturation should be avoided. 
The direct current servo motor with integral tachometer generator 
FD 84-21-1 built by the Diehl Manufacturing Company is considered 
representative of the drive motors available. The maximum speed of 
this motor is 415 radians per second and the maximum acceleration 

is 9,300 radians per second per second. For the selected gear ratio, 
the maximum controlled-member velocity is then 4,150 milliradians 
per second and the maximum acceleration 93,000 miliradians per 


second per second. 


Limitations on controlled- member velocity and acceleration 
establish a maximum error signal the systems investigated may receive 
as an input. In actual use of the systems no problem is presented. This 
is true since the system operates to keep the error signal nulled. In 
pulse testing, however, a limit is established on the pulse strength 
that may be used. Since pulse strength is a function both of pulse height 
and duration, it is possible to introduce too great a pulse height for 
the duration of the pulse. A pulse of this character will cause saturation 


of the drive motor velocity or acceleration or both. 


The relation between the motor inertia and the controlled-member 
inertia is an important parameter. In Chapter 3 it was shown that the 
controlled-member effective backlash is directly proportional to the 
effective inertia of-the motor and inversely proportional to the sum of 
the effective masses of the motor and the controlled-member. Witha 
gear ratio of one hundredth to one the effective mass of the motor when 
considered at the controlled-member is ten thousand times its actual 
value. In this study the effective mass of the motor is made equal to the 
effective mass of the controlled-member. This makes the effective 
backlash to be taken up by the controlled- member one half the backlash 
actually present inthe system. Also by this stipulation the actual inertia 
of the controlled-member is ten thousand times the actual inertia of the 
motor. It is evident that this inertia ratio is extreme but it is considered 
desirable. If the actual inertia ratio controlled-member to motor inertia 
is made smaller, the effective controlled-member backlash will also 


be reduced. To emphasize the effects of backlash the ratio of effective 
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inertias for this investigation will be taken as one half. 


Deviations in measurement of physical quantities will result. 
This class of deviations may be described as uncertainties because 
of their unpredictable nature. A maximum range of the output quantity 
may be specified which will include all the uncertainties. This range 
is defined as the uncertainty level. In the systems investigated, the 


uncertainty level is specified to be one half a milliradian. 
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СНАРТЕК 5 


PHYSICAL SYSTEM RESPONSES 


Physical feedback control systems using a motor and gear train to 
drive à controlled- member must exhibit non-linear performance char- 
acteristics for some operating conditions. For example, when the in- 
put quantities are large enough, the motor will exhibit either acceler- 
ation or velocity saturation or both, Backlash in the gear train will 
cause intermittant non-linear operation when the motion is reversed. 
Chapter 3 outlines a procedure for determining time responses of sys- 
tems with a little backlash. This method will now be applied to the 
physical systems described in Chapter 4. 

The investigation consists of three phases. The first is the cal- 
culation of ideal linear time responses for rectangular and displaced 
cosine pulse input functions. The second part of the study is the deter- 
mination of acttal responses from computed ideal responses when 
backlash is varied from one-fourth to one milliradian. The last phase 
is to evaluate the average error resulting from backlash in time re- 


sponses as well as frequency responses of physical systems. 
o-1 Assumptions 


The following assumptions are made to facilitate investigation. 
1. The average uncertainty due to interfering forces and tor- 
ques is considered to be zero. As will be shown later this 
uncertainty may be larger than the average error caused by 
backlash. The inclusion of the uncertainties due to interfer- 
ences in this situation obscures the error caused by backlash. 


2, The effects of friction are neglected, The controlled- 
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member velocity therefore remains constant during take-up 
of backlash. It is to be noted that Тизип '/ states that the ef- 
fects of friction are ameliorating. 

3. Bounce effects between driving and driven gears are not 
considered. 

4. Motor acceleration and velocity saturation will be pre- 
vented. To accomplish this, forcing function pulse strenghis 
will be limited. 

5. All initial conditions are zero. The motor gear is posi- 
tioned prior to test in contact with the controlled-member 
gear in the direction the pulse will drive the system when it is 
first applied. Controlled-member and mean displacement re- 
sponses will then be identical until separation first occurs, 

6. The amplitude of the controlled-member displacement 1s 
considered to be essentially zero when the mean motion dis- 
placement amplitude has settled to zero, This 1s a necessary 
assumption to establish convergence and integrability of the 
displacement response of the controlled-member. It is a rea- 
sonable assumption when the uncertainty level is comparable 


to the controlled-member effective backlash. 


272 Calculation of Ideal Linear Responses 





The Reeves Electronic Analogue Computer (REAC) of the Instru- 
mentation Laboratory, Massachusetts Institute of Technology, was 
used to calculate the ideal linear time responses, The instrumenta- 
tion of the REAC is described in Appendix B, Time responses were 
recorded on the output tables associated with the REAC. Output tables 
were used to obtain the accuracy required for graphical computation, 

Rectangular and displaced cosine pulse functions were applied to 
the AVSD and MAUSD systems, Pulse duration was one-tenth the nat- 
ural period of the AVSD system. Time responses were obtained of 
angular displacement, velocity and acceleration for various damping 
ratios, Angular velocity and acceleration responses were recorded 
for use in determining maximum pulse strengths for each system. 


These responses are shown in nondimensional form in Fig. 5-1 through 
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Table 5-1 Numerical values of first peaks of ideal linear 
response of AVSD and MAVSD systems to rectangular 
and displaced cosine pulse functions with 
(PD)(NP)R = 0.1 and PH = Imr. 


__ эт | тст 
AVSD MAVSD AVSD MAVSD | 
System system system system | 


Е | „2175 „2585 „1169 „1276 
И 0.7 . 2690 . 3030 . 1440 .1485 
EN 0. . 5020 . 5170 . 3065 . 2812 


ка 
o 


m 





‚0 21.65 21.65 13.15 
i 24. 25 24.50 14.12 
21 33.93 37.10 18.27 
‚ 0 3590 3092 
ДЫ | 3755 3205 
S 4325 3890 3875 
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Fig. 5-9. Numerical values from which these figures were consiructed 
are tabulated in Appendix С. Table 5-1 shows the numerical values of 
the response peaks for both systems. 

It may be seen from Table 5-1 that neither pulse function of this 
pulse duration is satisfactory as a forcing function, If a pulse strength 
sufficient to produce a ratio of one hundred to one peak displacement to 
effective backlash is used the resulting motor acceleration will exceed 
its limit and saturation will occur, It is therefore necessary to in- 
crease peak displacement response and investigate the effect on accel- 
eration and velocity peaks. 

Before discussing this part of the investigation, it 1s considered 
necessary to comment on the acceleration peak values caused by the 
rectangular pulse function. From Derivation Summary 3-1 it may be 
seen that the acceleration peak values should theoretically be indepen- 
dent of damping ratio and equal to и Ain: w^ for the AVSD system 
is 4007 radians per second per second. It may be seen from Table 5-1 
that the acceleration is dependent on damping ratio for the AVSD system and 
the value is not constant, Further it may be seen from Derivation 
summary 3-1 that the acceleration peak for the MAVSD system should 
be dependent on damping ratio and should vary from 1.1 we Ain to 
1.4 w^ Ain as damping is decreased from critical. w^ for the MAVSD 
system is 2800 rad/sec* . Referring again to Table 5-1, it may be 
seen that the MAVSD acceleration peaks for damping ratios of 1,0 and 
0.7 are approximately correct. The peak value for a damping ratio 
of 0.1 is not to be trusted at all. The cause of this bad reading is at- 
tributed to a scale factor omitted in the computer set-up or subse- 
quent calculations, It must also be pointed out that output table dy- 
namics could have contributed to these errors in computation since 
the acceleration peaks occur at time zero. 

From Derivation Summary 3-1 it may be seen that pulse dura- 
tion has no effect on the acceleration peak caused by a rectangular 
pulse function. This peak always occurs at the instant the system is 
disturbed. Its magnitude is directly proportional to the pulse height. 
It is evident therefore that no decrease in acceleration could result from 


increasing pulse duration if pulse height were kept constant, It may be 
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seen from Figs. 5-6 and 5-9 that the maximum acceleration occurs 
later in the response if a displaced cosine pulse is used. It was there- 
fore postulated that increased pulse duration with constant pulse height 
could produce decreased maximum acceleration if the forcing function 
were a displaced cosine pulse. 

The investigation was then continued using only the displaced co- 
sine pulse function as the forcing function, Only the seven-tenths 
critically damped form of each system is studied. Pulse duration is 
varied from one-tenth to one natural period of the AVSD system. The 
calculated responses in non-dimensional form are shown in Figs. 5-10 
through 5-15 for both systems. Appendix C tabulates the numerical 
values from which these figures were constructed. From these fig- 
ures it may be seen that for constant pulse height as the pulse dura- 
tion increases 

1. Angular displacement peak increases. 

2. The first angular velocity peak increases until the pulse 
duration is one half the natural period. "Thereafter it 
decreases. 

3. The second angular velocity peak increases. 

The first angular acceleration peak decreases. 

9. The second angular acceleration peak exceeds the first 
angular acceleration peak when the pulse duration is 
greater than three-tenths the natural period of the AVSD 
system. 

6. The second angular acceleration peak decreases after 
the pulse duration exceeds three-tenths the natural 
period. 

The investigation so far has shown that it is possible to increase ii. 
the controlled member angular displacement peak while decreasing the 
acceleration peak. 

Table 5-2 lists the numerical values of the response peaks of the 
AVSD system to a displaced cosine pulse function for various values of 
pulse duration, Table 5-3 lists the same information for the MAVSD 
system. From these tables it may be seen that the desired ratio of 
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Time - undamped natural period ratio T(NP)R 


Tig, 5-10 Displaced cosine pulse function input-angular chisplacement response 
for the AVSD system ‘for variation of the pulse duration (PD) 
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Table 5-3 Numerical values of peaks of ideal linear responses 
of the MAVSD system DR - 0.7 to a displaced 


cosine pulse function 
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angular displacement peak to effective backlash may be achieved. 
Pulse durations of seven-tenths the natural period of the AVSD system 
or greater may be used, Smaller ratios of displacement peak to ef- 
fective backlash may be obtained by using shorter pulse durations or 


smaller pulse strengths at the large values of pulse duration. 
023 Calculation of Actual Responses 


The second part of the investigation deals with the extrapolation 
to the actual response curves of the AVSD and MAVSD systems. This 
is accomplished using the method outlined in Chapter 3, The major 
portion of this phase is graphical integration and construction. In 
Appendix D will be found the complete calculation of a typical control- 
led-member angular displacement response curve, This section will 
consist only of a detailed description of the application of the method 
outlined in Chapter 3. 

The starting point is the ideal linear response peak table. Know- 
ing the acceleration and velocity limits of the motor and the gear ratio 
of the physical system to be tested, these limits are reflected through 
the gear ratio, These limits then become the maximum controlled- 
member velocity and acceleration permissable. These values for the 
AVSD and MAVSD system are respectively 4, 150 milliradians per sec- 
ond and 93, 000 milliradians per second per second, The next step is 
the determination of the maximum pulse height that may be used to 
disturb the system. In this investigation it can be seen that the accel- 
eration limit is the deciding limit, All calculations are made for a 
pulse duration equal to the natural period of the AVSD system. 

To determine the maximum pulse height, divide the allowable 
controlled- member velocity and accelerations by those tabulated in 
Tables 5-2 and 5-3 for the selected pulse duration. It may be seen that 
this maximum pulse height is 65.2 milliradians for the AVSD system 
and 64.5 milliradians for the MAVSD system. The mean member ang- 
ular displacement response may now be computed, This is accom- 
plished by dimensionalizing Fig. 5-10 for the AVSD system and Fig. 
9-13 for the MAVSD system is such a manner that the peak displace- 


ments are respectively 58.2 and 63.1 milliradians. The solid curve 


11 


in Fig. 5-16 represents the mean member motion of the MAVSD sys- 
tem, The next step is to dimensionalize Fig. 5-11 for the AVSD sys- 
tem and Fig. 5-14 for the MAVSD system. This is accomplished by 
making the mean member velocity first peaks 1524 milliradians per 
second and 1612 milliradians per second, The solid curve in Fig. 5-17 
represents the mean member velocity of the MAVSD system, The 
next step is to determine the controlled-member velocity curve, It 
will be recalled that friction effects were neglected in this investiga- 
tion. Therefore the controlled-member velocity is identical with 
mean member velocity until separation occurs, Separation occurs 
when the mean member velocity is a maximum ora minimum. At 
this instant of time the controlled-member velocity is at the maximum 
or minimum mean member velocity. The controlled-member contin- 
ues at this velocity until the area between the mean motion velocity | 
curve and the controlled-member velocity curve equals the controlled- 
member effective backlash, At this instant the controlled-member 
velocity assumes the mean member velocity. The controlled-member 
velocity curve for the MAVSD system with one quarter milliradian of 
effective backlash is shown in Fig. 5-17 as the solid curve including 
the shaded triangular areas. The hatched areas represent one quarter 
milliradian, The dashed curve of Fig. 5-16 is constructed by plotting 
the area of the hatched area of Fig. 5-17 as the ordinate difference 
from the solid curve of Fig. 5-16 from time of separation until the 
two velocity curves are identical. It is convenient to use a difference 
plot to show the variation of the difference ordinate. The difference 
plot for the MAVSD system is shown in Fig. 5-18 for an effective con- 
trolled- member backlash of one fourth millirad ian This difference is 
defined as (£)f (t). 

The method of determining the controlled-member angular dis- 
placement and velocity responses for one value of pulse height and 
one value of effective backlash for the MAVSD system has been ex- 
plained and illustrated, These steps must be repeated for each value 
of effective backlash at each pulse strength desired, No curves are 


included for other values of pulse height or effective backlash. 
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Fig. 5-16 Plot of angular displacement of mean member and controlled-member of 


MAVSD system for displaced cosine pulse function input showing effect of backlash. 
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574 Calculation of Average Percentage Error 





The final phase of this investigation consists of the calcualtion of 
the average percentage error in magnitude of the angular displacement 
of the controlled-member time response and frequency response, 


These functions are defined and derived in Derivation Summary 5-1. 


The angular displacement time response of a physical system 


with limited nonlinearities may be written in the form 


ADTTR ps) = f(t) + (E)f (t) | (1) 


where ADTTR angular displacement transient time response 


(ps) 
of the physical system 


f(t) = angular displacement transient time response 
of the ideal linear system with same paramet- | 
ers as the physical system but with no non- 
linearities 

(€) f(t) = angular displacement deviation of ADTTR 


(ps) 
from f (t) 


Each quantity defined above is a function of time (t). (£)f (t) is calcu- 
lated as outlined in Chapter 3. 
(а) Angular displacement time response 


(9) 


Cramer shows that Schwarz inequality may be written in the 


form 


"e p 
| | cerco "ass 4 2+2 Keoee [ee] (ежој | dt (2) 


со 
Neglect the term | |©ге)| dt as small compared to other terms, 
-00 . 
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Derivation Summary 5-1. Derivation of error describingfunctions for 


a physical system with limited non-linearities. (Page 1 of 5.) 
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00 2 
Subtract | |“) dt from both sides of equation (2) 
-00 


| [EOE EE] at- Fa) ac |2 («ӘЛІ КӘ» (3) 
-00 -00 ! -00 


oo 
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= 00 = 00 -00 








« 
со oo 
| ік) ? gt f [eeo] dt 
-00 00 
(4) 
А со 00 
Ше 100 { есеј dt 100 ( еса 
-00 =00 
APEM tr) > (5) 
со 
| го а 
-00 
where 
APEM ttr) “ average percentage error in amplitude of the tran- 
^. sient time response of a physical system 
Therefore l 
со 
200 || | (ec || teo | at 
-00 
APEM in percent (6) 





< 
tp) 

f(t) dt 
мо 


(b) Expression of average percentage error in amplitude 


of the transient time response of a physical system 


Derivation Summary 5-1 Derivation of error describirg functiors for 


a physical system with limited non-linearities. (Page 2 of 5.) 
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' 


The angular displacement frequency response of a physical SyS- 


tem with limited non-linearities may be writteniin the form 


- 1) 
ADFR(, * F (8) + (€) F (8) (7) 
where 
aE (ps) = angular displacement frequency response of 
the physical system 
F (8) = angular displacement frequency response of 


the ideal linear system with same parameters 
as the physical system but with no non-linear- 
ities 

&) к (В) = deviation of ADF'R ps) from F (8) 





В = frequency ratio and is the forcing frequency 


(c) Angular displacement frequency response 


divided by some arbitrary reference frequenc 
Applying the steps outlined in section (b) it can be shown that | 


100 | | FHE O|" aB - 100 || | (8)| ав гоо\ (ЕР (8) |е(в)јав 
= 00 = 00 = 00 


< 





00 2 m oO 
f |= (812 ав f |к(8]° ав 
m -00 
(8) 
The mathematical expression F(8) is used to describe a physical re- 
lationship. In nearly all physical systems F(f8) is an even function, 


Negative frequencies do not exist. Therefore 


Derivation Summary 5-1 Derivation of error describing functions for 


a physical system with limited non-linearities. (Page 3 of 5,) 
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100 ( | F(8-€)F (8) 46-100 ( | (8) ав 400 | Kr (8)| [F (6)]8 

– 00 -00 0 

2 —— D 

од 2 со 2 
Е à F 

\ [rotas М (8)|“ ав 











5 oO oO 
iine 100 f |F EEF E| 48-100 | [POl ag 
1 ~ со v —00 





4 


r [ғ (8)[ ав 








| where 
APEM fr) = average percentage error in amplitude of the fre- 
quency response of a physical system 
QO 
К оге 200 ( |ЕЈЕ(В)| |Е(В)|а8 
АРЕМ(,;5 in percent (11) 


f [rey ав 
0 


Ae ROO ^ ти ер ФЕР шры: f 


(d) Expression of average percentage error in amplitude of 
the frequency response of a physical system 


Morse and ок 


ы 


show that Parseval's formula may be | 
written 


| Fy] Lo С е? dt (12) 
-00 -00 


From Equation (12) it may be seen that 


Derivation Summary 5-1 Derivation of error describing functions for 


a physical system with limited non-linearities. (Page 4 of 5.) 
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o0 2 со 2 со 2 00 2 
f [peer [Гав - | [F (6) ав f (+ (| а‹ - ( е [а 
-00 -00 -00 =00 





Ü |0812 ав С eo" at 
Е -00 


Therefore 


АРЕМ, 





- АРЕМ 


fr) (ttr) 


(e) Relationship between average percentage error in amp- 
litude of the transient time response and average per- 


centage error in amplitude of the frequency response of 


a physical system 


Derivation Summary 5-1 Derivation of error describing functions for 


a physical system with limited non-linearities. (Page 5 of 5.) 


The deviation plots described in 5.3 are plots of (É)f(t) as stated. 
Plots similar to the solid curve of Fig. 5-10 are constructed for each 
value of pulse strength to be considered for each system. A plot of 
2| (£)£(t) ||f(t)|is then constructed for each effective backlash value by 
multiplying the ordinates of the (€)f(t) and f(t) curves together. A plot 
of «|^ is constructed for each pulse height for each system. These 
curves are integrated, Graphical methods were used by the authors, 
By combining the integrated functions as shown in section of Deriva- 
tion Summary 5-1 APEM ttr) is obtained. 


Fig. 5-19 is a nondimensional plot of APEM as a function of 


(ttr) 
pulse strength and effective backlash for the AVSD system, Fig. 5-20 
is a similar plot for the MAVSD system. The quantities used to non- 
dimensionalize these curves are effective backlash ratio- BR- and ref- 
erence pulse strength -( PS) ref. These two quantities are defined be- 
low, 
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The curves of this 
fia. are for 
EISE (РОММР) =1.0 (DR)= 0.7 
L1 A1 ]- 1-2] E 
LI LA] TO ш 
ПГ 1 
а ae 
ЕРРИЕЕЕМШЕГ | 7777772 
а Т ые ШЕШ d 
ој РК У | | Ра РАК 7" 
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coc LIN Ne a tt 
LLL EAM L NL HEEL LLL LL LL 
HH 
А). том 


APEM (itr) 





0.4 0.6 0.8 І.0 


PS 
к Strength Ratio (9S) e 


Average percentage error in transient time response 


Fig. 5-19 Angular displacement average error produced by backlash 
in the AVSD system with DR=0.7 and (PD)(NP)R = 1.0 
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Average Percentage Error in Transient Time Response APEM atto) 






jt tt (Note: The curves of this fig.| 
ШЕШ “ЖЕ are for P 
ШЕШ _ || (РрХМРЈЕ. + 1,0 (DR)=0.7 | | 
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ДО ДЕ ве ШЕШ 
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TORAN 
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ENEN 
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А0 
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Pulse Strength Ratio (PSR) = 





E 
Fig. 5-20  Angular displacement average error produced 


by backlash in the MAVSD system with (DR) 
0.7 and (PD)(NP)R s 1.0. 
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д liom) n 


BR - 
А Т 
т п 
where 
Ы (ст) controlled-member effective backlash 
n = gear ratio motor to controlled- member 
• 
An - maximum motor velocity 
T, - natural period of AVSD system 
(PD) (PH) 
(PS) - | осе (тсе eme dt 
ref 2 о) 
0 п 
мћеге (РН) = pulse duration 


(РН) max ^ maximum pulse height possible avoiding motor 


saturation 
t = time 
Sa natural angular frequency of AVSD system, 


Table 5-4 lists the numerical values from which Fig. 5-19 was con- 
structed. Table 5-5 tabulates the same information for Fig. 5-20. 


From Derivation Summary 5-1 it may be seen that АРЕМ у) 


is equal to APEM e. . Therefore although not so labelled, Figs. 5-19 


and 5-20 show the variation in АРЕМ у) for different values of BR 


and (PS) (PS) for the AVSD and MAVSD systems respectively. 
ref 


The effect of the uncertainty level in a system is handled in an 
analogous manner. It is apparant that the uncertainty level is present 
continuously throughout the motion. For this case, (&)f(t) remains 
constant and is equal to the magnitude of the uncertainty level (U)f(t). 
To compare the effect of uncertainties and backlash, a sample com- 
putation of the average percent error due to an uncertainty level of 


one half milliradian in the MAVSD system is included in Appendix D. 
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Table 5-4 Average percentage error in transient time response 
of the AVSD system due to backlash. 


| 10 | 07 | os | ozs | 


APEM APEM APEM APEM 











PSR 


А(ы) (ст) 
(mr) 










Table 5-5 Average percentage error in transient time response 
of the MAVSD system due to backlash. 





| no | oz | os | os — 


APEM APEM APEM APEM 
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The results indicate that an average error of 2.42 percent will exist 
due to the uncertainty level. Fig, 5-20 shows the average error for a 
total backlash of one half milliradian to be 0.798 percent for a pulse 
strength ratio of 1.0. It is concluded from this comparison that when 
the magnitudes of the backlash and uncertainties are equal, the re- 
sponse deviations due to backlash are negligible. 

Investigation of the technique of pulse distrubance testing of feed- 
back control systems discolses several significant considerations. 
Velocity and acceleration limitations of a given system must not be 
exceeded, The peak amplitude of the displacement response of the 
system must be large relative to the backlash. These requirements 
can be met simultaneously by proper choice of the pulse function. In 
this connection, pulse shape and duration, as well as pulse strength, 
are important. Utilization of the correct pulse function minimizes 
the response deviations due to backlash, and these deviations are neg- 


ligible for nominal amounts of backlash. 


J 
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СНАРТЕК 6 
CONCLUSIONS AND RECOMMENDATIONS FOR FURTHER STUDY 


Responses of physical systems to pulse functions are used for 
evaluating the linear properties of the systems. The intent of this 
study is to determine how to apply pulse function inputs which yield 
information about the linear properties of systems having some signi- 


ficant non-linearities. 


Physical systems may exhibit a wide range and variety of non- 
linearities. In order to develop methods for studying such systems, 
the non-linearities have been restricted to backlash in the gear train 
and saturation of the torque source. Velocity saturation exists when the 
torque source output member has a maximum velocity it cannot exceed. 
Similarly, acceleration saturation implies the existence of an upper 
acceleration limit. These are quite general since every physical system 
possesses maximum velocity and acceleration limits. The procedures 
described in the previous chapters are adaptable to other types of non- 


linearities. = 


Two representative physical systems have been studied nere. 
The AVSD system ‘is a second order type, characterized by angular 
velocity signal feedback damping. Modifying the angular velocity feed- 
back signal converts it into the third order, MAVSD system. Non- 
linearities in the form of small amounts of bac klash as well as velocity 
and acceleration limits on the motor are introduced in the ideal systems. 
Rectangular and displaced cosine pulse function inputs are applied and the 
system responses studied. The effects of pulse shape, duration and strength 
on the system responses are examined. The following conclusions are 


derived from this study. 
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6-1 Conclusions 


(a) Velocity and acceleration limits of the motor must not be ex- 
ceeded. Saturation of any component can disto rt the transient response 
beyond usable limits. Correlation with the ideal responses determined 


from linear analysis is then impossible. 


(b) The pulse function input must have continuous derivatives to 
avoid infinite inputs. This is first described in the theoretical analysis 
of Derivation Summary 3-1 of Chapter 3. The components of a physical 
system are incapable of accepting an infinite input. They may be driven 


into saturation by such an input. 


(c) The pulse strength must be adequate to force an output large 
compared to backlash and uncertainties. A ratio of the peak amplitude | 
of the transient response to the effective controlled- member backlash 
of one hundred to one is desirable. If this is achieved without saturation 


distortion, the deviations from the ideal response are negligible. 


(d) The pulse duration should be approximately equal to the undamped 
natural period of the dominant mode of the system. Under this conditi on 


the pulse strength may be adequate without inducing saturation. 


(e) A rectangular pulse function is unsuitable for the systems 
in this study. This is discussed from a theoretical point of view in 
Derivation Summary 3-1 of Chapter3 and verified numerically in Chapter 5. 
The discontinuous nature of this function and its derivatives makes it 
impossible to avoid saturation of the acceleration and velocity for the in- 


vestigated systems. This condition is true regardless of pulse duration. 


(f) Functions such as the displaced cosine pulse may meet the 
requirements stated above. They possess continuous derivatives. Further, 
it is possible to introduce adequate pulse strength by adjusting the pulse 


duration without inducing saturation. 


(g) The average non-dimensional error in amplitude due to backlash 
and uncertainties is the same in both the time and frequency domains. This 
is shown in Derivation Summary 5-1 of Chapter 5. Since the correct 


choice of a pulse function results in a negligible average amplitude error 
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in the transient response, the average error in the amplitude of the 


calculated frequency response will also be negligible. 


(h) When the backlash and the uncertainties are of the same order 
of magnitude, the average amplitude error due to backlash is smaller 


than that due to the uncertainties. 


6-2 Recommendations for further study 





The investigation reported in this thesis is limited to a study of the 
non-dimensional error in the frequency response amplitude. The effect 
of backlash, uncertainties and saturation upon the dynamic response 
angle has not been determined. However, the necessary information 
can be obtained by an extended study of the kind used herein. It is 


recommended that this be done. 


Further study of the problem of calculating the frequency response 
from the transient response is suggested. Cognizance of the character 
of transient response errors arising from backlash and other uncertainties 
is important here. The distribution of errors in the frequency response 
may be ascertained from carrying out this recommendation. 

A theoretical investigation of the sinusoidal responses of systems 
having backlash and saturation effects will be useful in comparing direct 


frequency response testing methods with pulse techniques. 


~ 


ж 


Experimental verification of the conclusions described above is 
strongly recommended. More convincing conclusions may then be drawn 
as to the relative merits of pulse techniques versus direct frequency 


response testing methods. 
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APPENDIX A 


DERIVATION OF THE PERFORMANCE EQUATIONS OF THE SYSTEMS 
INVESTIGATED 


The performance equations of the positional servomechanisms 
described in Chapter 4 are derived by summing the torques acting on 


the controlled-member. 


A-1.Positional Servomechanism with Angular Velocity Feedback 


Damping -AVSD System 


Consider the components of the functional diagram of Fig. 4-1. No- 


tation used is defined there. 


(a.) (Sc) = angle signal comparator 


"(sc) З(зс)[ГА;е im) j А (сто) 


*(sc) ^ sc) [A;e ) З(сто) (A-1) 
(b.) (pa) - кыч рет 
(ра) > (ра) [е;е | °(зс) о 
(c.) (vsc) Velocity signal comparator 
"(vsc) | (vsc) Ге; е] [ (ра) | е (ува) а) 
(d.) (da) = Drive amplifier 
(А-4) 


*(da) ^ Э(аа) [е;е] *(vsc) 
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(е.) 
(tg) = Torque Generator 


Мағ) 7 igt) Pen: Ars] (tg) fe; M] *(da) 


: E. 
^? (80 [Aca Atsm| За) А; м] Atem) (А-5) 


Reduction of the output of the torque generator indicated in equation (A-5) 


is the result of motion of the torque generator rotor relative to its case. 


(f. ) 
| (avsg) = Angular Velocity Signal Generator 


"(avsg) 5 S(avsg) | A;e P (tsm) 


_ 2 | s | 
"(avs)" "teo [Aou Atsml З(амзе) |Аге] ^em) — 0^9 
(g.) 
(уза) = Velocity Signal Amplifier 
(ува) | S(vsa) (еҙе| "(avsg) x) 


Combining the above relationships for the inner and outer loops of the 
functional diagram, Fig. 4-1. yields: 


NC F Slet Aem Atam Уа) [е; М] *(da)[e;e] Э(узс) [е;6] (ра) [е;е] 
2 : 
(sc) fa; Alem) - Sigt)  omiAtem} lavse) [Arq S(vsa) [е;еј 


А 2 5 : 
"tvse)[e:;e] "(da)[e;e] Sit g) ем А (ст) - Teo [^ cA su] ejoa] Қ (ста) 
(А-8) 


92 


An inertia reaction torque is present when the controlled- member is 


accelerated with respect to inertial space. 


Мат) * 7 Iemyteo) ^ft - em] (А-9) 


If the base does not move with respect to inertial space, 


^r = cm) ы Ж” 


[ES А ВВ 

Gr) 7” Қош)ев) ^ (em) n 
A basic dampingiorque exists due to the relative motion between the con- 
trolled-member and the supporting base. It is partially dueto friction 


in the bearings and gear train, and to air damping. 


Ма) =- С TS (8-11) 
This reduces to 
Ма) = - Са P (em) (А-12) 


ж 


if the base does not move with respect to inertial space. 


~, 
a 


= 


Torque Summing Equation at the controlled- member 


Mag) * Маг) 7 Миа) 7 Ману) 7 0 (А-13) 


Substituting equation (A-8), (4-9), and (4-12) in equation (A-13) gives 
the following performance equation of the AVSD system: 


Sí gt) КОШ tem] P(tg) le; M] (da) [e;e] (увс) [е;е] S( pa) Бље] (вс) [А; 4 (О) P Cem) 


2 i . Е 
F [A cni tem! Лого ГА еј “беде еј З(увс)Ге;6] (да) (ese; Ate) le: М] 


B ЗЕ "lu e А 


| 


А(ста) + М = 0 (А-14) 


cemo а (intfr.) 
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Define: 


S(ps) A;M ` "G9 [^en] 5+а)[е; М] Э(аа)[е; е] (уве)[е;е] (ра)[е;е] 


elastic coefficient (A-15) 


"(sc)[A;e] У 


S(rb)(av)|A; M] = Э(аузв)[А; < (уза)[е; | (vso) [e:3 ^ta) [ed "te)[e; M] 
(A-16) 


; Е 2 | 
“(es) [asm | (et) (Аста “tem ЕСЕ мј" Қалауы) Е 
(А -17) 


Substitution of equations (A-15) through (A-17) into equation (A-14) 


and rewriting the result in operator form yields: 


2 : = 
[emiten P * Sos)[A; M] 1 Alem) ^ Stps)[A; Mf ^ (em) * Маім) 
(A-18) 


A-2. Positional Servomechanism with Modified Angular Veolcity Feed- 
back Damping - MAVSD System. 


Consider the components of the functional diagram of Figure 4-2. 
Note that the only difference in components from those of Figure 4-1 
is the addition of a velocity signal modifier in the angular velocity feed- 
back path. The velocity signal modifier is a simple passive network 


differentiating circuit which has the following performance function: 
(СТ) ysm)? 


It (CT) тур 


FF (vsm)fe;e] . (А-19) 


With the velocity signal modifier in the angular velocity feedback path 


the torque generator equation is as follows: 
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(СТ) өт)? 2 


Mtg) З(рв)[А;М] С) (ст) 7 1+ (СТ тур | "(89 [cs cse] 


"^ CH Ж 
*(to)(av)[A:M] “(cm) ~ БОГА А] ЕМ] У (сто) 
(А-20) 


The inertia reaction and damping torques derived in equations (A-10) 
and (A-12) are unaffected by the velocity signa! modifier. Subsituting 
the expressions for the torques into the torque summing equation (A-13), 
simplifying and writing in operator form produces the performance 


equation of the MASVD system. 


I MOM. TO vem? d 
(emea)? * Pte)[A, iA LL] Pg [A;M] 11 (Ст) ав 19 v[A; M] P Alem) 
7 5 vsm 


Зар АМ A (em) * (inttr) (A-21) 


Attention is invited to the fact that the coefficients of equations (A-18) 

and (A-21) may be kept constant while permitting the sensitivity of the 
gear train to vary, This is accomplished by adjusting amplifier sensi- 
tivities where necessary to compensate for a change in gear train 
sensitivity. As a consequence, it is possible to use the gear train 
sensitivity as a variable parameter in the study of pulse testing techniques 


for systems with backlash. 


A-3 Numerical Performance Equations for the AVSD and MAVSD Systems 


Numerical values for the coefficients of the performance equations 
derived above were selected from Lees 9) for a system damping ratio 
of seven-tenths critical at high frequencies. Coefficients were then 
determined for critical damping and one-tenth critical damping without 


changing the undamped natural frequency. Interference torques were not 
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considered in this investigation. These equations are shown below. 


AVSD System - undamped natural frequency w = 63.3 rad/sec 


n(ps) 
(1) Critical Damping, (DE) 5s) - ]1.0 
(p^ + 126.6p + 4007) A = 4007 m in) (A-22) 


(cm) 


(2) Seven-tenths critical damping, (PR) ps) 7 0.7 


(p^ + 88.7p + 4007) Alem) = 4007 B (4n) (A-23) 
(3) One-tenth critical damping (DR) 55) - 0.1 
2 2 
(р + 12.66р + 4007) Aem) - 4007 А n) (A -24) 


MA VSD System - undamped natural frequency w - 53.6 rad/sec 


n(ps) 
(СТ) sm) - 0.0993 sec. 


(1) Critical Damping (DE) 5s) - 1.0 
3 2 
t 127.4p + 4030р + 40350) А 40350 + 
(р р р ео =. 


(А-25) 
(2) Seven-tenths critical damping (БА) ба) - 0.7 


3 2 
(p + 99.4р + 4030р + 40350)А = (40350 + 4010р)А 


(іп) 
(A -26) 


(cm) 


(3) One-tenth critical damping UR) S) ~ 0.1 


2 
(р? + 22.83р + 4030р + 40350)А = (40350 + 4010р)А 


(іп) 
(А -27) 


(ст) 
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APPENDIX B 


REAC INSTRUMENTATION 


The instrumentation of the second order linear differential equa- 
tion of motion of the AVSD system for REAC simulation was accom- 
plished using standard techniques. For this reason, a discussion of 
the problem is not deemed necessary. The third order, linear dif- 
ferential equation for the MAVSD system required a somewhat differ - 
ent approach to effect correct REAC simulation due to the presence of 
the derivative of the input quantity. The method of instrumentation is 
developed below and a diagram of the basic REAC setup is shown in 
Fig. B-1, 

The MAVSD performance equation, as derived in Appendix A, is 


of the following general form 


3 2 
(р +a,p +a,P * ag) А (ст) = (bo tbi p) A Gn) (B-1) 


To obtain the solution of equation (B-1), where Ain) represents the 
desired pulse function input, instrument the following equation by 
standard REAC methods. (See Fig. B-1) 


3 2 i = ~ 
(р фар +ајр + aq) A oe hg Adn) (B-2) 


Let 
À ün) 
o E е. (B-3) 


3 2 
(p а р +а,р+а‹у) 
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"urejsÁs ASAVNWN 993 toy dmas DVWH oTSegH Т-Я ‘814 


"uroJ0op stu u1 pepniou| 100 зилоб 2әлуцдшо рио 20042 21095 әш : 2403 


HD о иу за > | ien 


шә tan 
AD 
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БК aum ваз 





атаў! 
104100 шоу айт 






31871 31991 
LNALNO шоу от (айо шә. д 0p 


Then from (B-1), 


А (ст) = by y +b, py 


(B-4) 


The first and second derivatives of equation (B-4) in operator form 


are. 


Р Ы) 


2 2 3 
pg У ал о 


From equation (B-2), 


= 


А tem) 


3 2 
(p та P +a, Pt ao) 


bg A 


The derivatives of (B-6) are, 


p А tem) 


2 1 
DES (cm) 


Then, 


Substituting in equations (B-4) and (B-5) gives: 
A = А! Шы 
(стп) (cm) ' bg 


2 ! 
pa = РА 


bg P Y 


2 
bg P Y 


2 
bgPYtbP y 


b 
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(in) 


р А! 


(а) (В-5) 


(b) (B-5) 


(B-6) 


(а) (8-7) 


(б) (В-7) 


(а) (В-8) 


(b) (B-8) 


(a) (B-9) 


(b) (B-9) 


b 


1 30, E 
(cm) * bo Р А (ст) ӨЛЕДІ 


2 к. 
р А (сто) р А 


Multiply the equations of (B-9) by the appropriate coefficients from the 


characteristic equation of (B-1). 


= ap A! (a) (B-10) 


Ш 1 
80^ (cm) (cm) * 0 bg pA (cm) 


а | РА (сл) = а; рае) +a p P (b) (B-10) 


2 2 bi з 
' (B-10 
ар А а, р A! (cm) * 325: b. р А (em (c) (B-10) 


(cm) | 
Note that in the basic REAC setup of Fig. B-1 the quantities (ay А! (стају 


апа (аур е) are present. Тћеп to determine (ао А ) from 


(cm) 


equation (B-10)(a), itis necessary to multiply (а |) by a factor (ki) 





such that, 
by 
a Kk, = а 
1 1 0 50 
ог agb, 
ky = E (a) (B-11) 
Similarly for (B-10)(a) and (b), 
een 
_ \ 4 
ko = Е (b) (B-11) 
270 
aab 
ч (c) (B-11) 
3 bo 


By the above procedure of instrumenting equation (B-2) and cor- 
rectly summing quantities at pertinent points, it is possible to obtain 
the displacement, velocity and acceleration responses of equation 
(B-1) to a desired pulse input to the system, The method avoids 
REAC differentiation of the pulse input and gives reproduceable time 
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function responses. 

Rectangular pulse generation was accomplished using a step in- 
put which was returned to zero magnitude at the end of the desired 
pulse duration by a relay. Displaced cosine pulse function generation 
was simply obtained by intrumenting an undamped second order sys- 
tem, Pulse amplitude was controlled by varying the initial condition 
on the second integrator. Pulse duration was established by a relay, 
and the undamped natural frequency was adjusted through proper 
choice of equation coefficients to give one cycle of the oscillation for 
the period of the pulse duration. Scale factors to meet loop gain lim- 
itations of the REAC were such as to make any minor deviation in 


pulse shape insignificant in their effect on system response. 
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APPENDIX C 
TABULATION OF NUMERICAL VALUES OF T HE RESPONSES FOR 
IDEAL LINEAR AVSD AND MAVSD SYSTEMS 


Table C-1 Values of controlled-member displacement ratio for rect- 
angular pulse function input to AVSD system (PD)(NP)R = 0.1 





T(NP)R 

0.000 0.000 0.000 0.000 
0.050 0.095 0.090 0.079 
0.151 0.654 0.462 0.390 
0.252 0.978 0.523 0.420 
0.352 0.915 0.420 0.333 
0.453 0.524 0.267 0.235 
0.554 0.000 0.135 0.155 
0.655 -0.472 0.045 0.099 
0.755 —. -0.714 0.000 0.061 
0.856 -0.671 -0.021 0.037 
0.957 -0.388 -0.024 0.022 
1.057 0.000 -0.019 0.014 
1.158 0.363 -0. 011 0.009 
1.260 0.518 -0.006 0.006 
1.360 0.488 -0.002 0.004 
1.460 0.287 0.000 0.003 
1.510 0.157 0.000 0.002 
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Table C-2 Values of controlled- member velocity ratio for цестанеш И 
(PD)(NP)R = 


pulse function input to AVSD sysiem 






T(NP)R 


каше Se See ее ее 
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ref - ^em(max)(DR = 0.1) 


0.1 


\ 





.000 
.430 
566 
213 
‚081 
. 149 
.135 
.103 
Um 
‚046 
.029 
.016 
.013 
.008 
„006 
.005 
.002 
.001 


Table C-3 Value of controlled- member acceleration ratio for rectan- 


gular pulse function input to AVSD system (PD)(NP)R = 0.1 


Sk AT: CA OR (EA ПЕЊЕ РАЈ 
| 


T(NP)R 


кі юе кі ра e eae о O СО O O со о 


‚000 


025 


„095 
121 


201 
302 
403 
903 
605 
705 
806 


.907 
.007 
.108 
208 Т. 
310 
‚410 
‚510 


Быны саннын да 





.314 0.000 0. 
‚051 0.057 0. 
‚208 0.057 0. 
.357 0.051 0. 
.359 0.033 0. 
235 0.020 0. 
.038 0.006 0. 
.152 0.000 0. 
‚258 0.000 0. 
‚268 0.000 0. 
177 0.000 0. 
‚033 0.000 0. 
Aref | Bom (max)(DR = 0.1) 
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Table C-4 Values of controlled- member displacement ratio for rec- 
tangular pulse function input to MAVSD system (PD)(NP)R = 0.1 


T(NP)R em ret 
DR = 0.1 
“а 
0.000 0.000 
0.063 | 0.152 
DOE d -- - 
0.126 | 0.514 
0.189 | 0.817 
0.246 т en. 
0.252 0.976 
0.258 | шақы. 
0.289 | 1.000 
0.314 ГЕ о: 
0.315 о ----- 
0.327 | 0.967 
0.378 | 0.837 
0.503 | 0.270 
0.630 | -0.339 
0.693 -0.556 
0.755 -0.664 
0.794 | -0.679 
0.882 -0.560 
1.007 ОНО 
1.132 0.263 
1.126 0.462 
1.285 0.490 
1.385 0.405 
1.510 0.118 
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ref ^ Зет (max)(D R = 0.1) 
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Table C-6 Values of controlled-member acceleration ratio for rec- 








tangular pulse function input to MAVSD system (PD)(NP)R = 0.1 
T(NP)R | A em! Pret 
DR = 0.1 | DR = 0.7 DR = 1.0 
i 
0.000 0.000 0.000 0.000 
0.013 0.405 0.127 № № = 
0.025 1.000 | ---- 0.200 
0.045 1 ----- 0.279 (|| --- 
0.063 scene 0.206 0.122 
0.076 | ----- 0.110 “fl --Е 
0.101 0.716 -0.022 0.064 
0.113 | ----- -0,178 PE u E 
0.126 =0.378 > (| == И По и 
ooa a -0.163 
0.189 | -0.505 -0.106 -0.092 
0.252 | -0.550 -0.063 -0.047 
0314 №. -0.507 № Ш 8 ---====_ (а 
0.378 | -0.394 -0.014 -0.006 
0.503 ` | -0.047 0.007 0.004 
0.630 | 0.270 0.012 007 
0.755 | 0.380 сз 006 
0.882 | 0.270 0.007 — NP S 
1.007 | 027 | EE 0.003 
| 1.132 | -0.195 0.007 | ----- 
| 1.260 | -0.270 ---- 0.001 
1.385 | -0.195 |2 ----- | ----- 
1.510 | -0. 014 -0.001 | 0.000 


ео НИ 





ref - А ср (max)(DR - 0.1) 
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APPENDIX D 


TYPICAL CALCULATION OF ERROR 


The purpose of the following computation is to illustrate the ap- 
plication of the methods of Chapters 3 and 5 to determine the average 
response error resulting from backlash. Tables 5-4 and 5-5 sum- 


marize the results of similar calculations. Figures 5-19 and 5-20 


show the variation of the average percentage error with pulse strength. 


m MAVSD System 
1, Pulse Function: Displaced Cosine 


(PD)(NP)R = 1.0 


2. Motor Limitations: At the torque summing member,maximum 


acceleration = 9300 rad/, c? 


‚ж 


maximum velocity - 415.0 rad/sec 


3. ~“ Traain: P (tsm) /; - 100 
в (cm) 
Total backlash = 0.5 milliradian = Ау) 
Effective backlash = 0.25 milliradian "ASIDE 


4. Pulse Strength: 
(a) From Table 5-3, a pulse height of 1 mr. gives 


А (ст) = 1441 MP / ес“ 


(b) From motor characteristics and gear ratio, 
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.. 





ах) 93000 таг. / 2 
93000 . 
(c) ЫН)... =a earn 64.5 Mr. 
PH)(T_) 
(d) (PS) = 5 H Е = 3.21 mr.sec, 


0, Peak amplitude of displacement response: 
From Table 5-3, 


E Cm) ende) - (PH)x 0.967 


ах) = 64,5 x 0,967 = 62,4 mr, 


6. Peak amplitude of velocity response 


From Table 5-3, 


A) - (PH) x 25 = 64.5 x 25 = 1612 mr, /sec. 


From motor characteristics and gear ratio, 


сий л) 4150 mr sec. 


Therefore the velocity is well below it's limit, 


7. Application of the effective backlash to the velocity response 


of the mean point, (Shaded triangular areas of Fig. 5-17.) 


5 AVel x Time 
2 


А (ы1) (ста) 2280259 mm 


Note: As the velocity peaks decrease in amplitude, a great- 
er period of time is required to take up the effective 
backlash, 
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8. Deviation due to backlash Fig. 5-16 


(a) During time the backlash is taken up the deviation rises 
parabolically. 

(b) Deviation remains constant at the value of the effective 
backlash until the velocity reaches a negative maximum. 

(c) Deviation falls parabolically to zero as backlash is 


taken up in the opposite direction. 


9. Average error in percent 


oO 
2 | | ecol H1 at 


-00 
АРЕМ © EL D 
o0 
f еу dt 
-c0 
where (£)f(t) = ordinate of deviation curve, Fig, 5-18 
f(t) = ordinate of displacement curve, Fig. 5-16 


(a) Integrating graphically, 


f КОЈЕ = 173 (mr)? sec 


-00 


р со 2 
i 2 | ОП (Б а = 1,38 (mr)“sec 


-00 
1, 38 
173 


x 100 = 0.798% 





(b) APEM Цу) 


10. Average error for uncertainty level (U)f(t) of 0, 5 mr. 


со 
2 ( КӘ себе! = 4.18 (mr)“sec 
c =O 
„а аса 
APEM r) = 7179 Х100. = 2,42% 


GLOSSARY OF SYMBOLS AND NOTATIONS 


Self-defining notation developed by Draper et a1) has been used 
throughout this thesis. 


ADFR - 


ADTTR - 


‚ж АРЕМ E 


AVSD - 


bl(ecm) - 


amplitude 

angular displacement 

angular velocity 

angular acceleration 

angular displacement frequency response 
angular displacement transient time response 
average percentage error in magnitude 
angular velocity signal damping 
non-dimensional frequency ratio 
backlash 

effective controlled- member backlash 
backlash ratio 

correction 


controlled- member 
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DAR 
DCPF 
DRA 


DR 


eq 
(€ )£(t) 
F( 2) 
f(t) 

fr 


ER 


id 


IFF 


MAVSD 


M 


dynamic amplitude ratio 

displaced cosine pulse function 

dynamic response angle 

damping ratio 

natural logarithm base 

equivalent 

error in the function of time 

working variable form for frequency response 
working variable form for transient time response 
frequency response 

frequency ratio 

inertia 

identified 

input forcing function 

imaginary component 

modified angular velocity signal damping 


torque 
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п 
ге! 


OR 


PF 


PD 


(PD)(NP)R 


ps 


PSR 


REAC 


milliradian 


frequency in cycles per second 


forcing frequency in cycles per second 


reference frequency in cycles per second 


output response 


derivative operator 


performance function 


pulse duration 


pulse duration natural period ratio 


pulse height 


pulse strength 


positional servomechanism or physical system 


pulse strength ratio 


quantity 


Reeves Electronic Analogue Computer 


sensitivity - dimensional similarity factor 


time 
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Т 
n 
(T)(NP)R 
an 
tsm 
ttr 


(0:06) 


We 


гө. B 


undamped natural period 

time natural period ratio 

working variable expression for (PD)(NP)R 
torque summing member 

transient time response 

uncertainty in the transient time response 
angular forcing frequency 

angular undamped natural frequency 
working variable form for DAR 

working variable form for DRA 


working variable form for 5 
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